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We have experimentally studied the non-equilibrium transport in low-density clean 2D electron 
systems at mesoscopic length scales. At zero magnetic field (B), a double-peak structure in the 
non-linear conductance was observed close to the Fermi energy in the localized regime. From the 
behavior of these peaks at non-zero B, we could associate them to the opposite spin states of the 
system, indicating a spontaneous spin polarization at _B = 0. Detailed temperature and disorder 
dependence of the structure shows that such a splitting is a ground state property of the low-density 
2D systems. 
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Spin polarization of electrons in low dimensional sys- 
tems at S = has recently attracted extensive theoret- 
ical and experimental attention. In the metallic regime, 
this is commonly attributed to the spin-orbit (SO) ef- 
fect, where the bulk inversion asymmetry (the term) 
and the interface inversion asymmetry (the Rasba term) 
result in an energy separation of several millivolts be- 
tween the spin bands in narrow band semiconductors Q . 
At low electron densities however, the electron-electron 
interaction dominates over the kinetic energy, and a ho- 
mogeneous 2D system becomes unstable to spontaneous 
spin polarization (SSP) due to exchange 2, 3]. Since the 
influence of interaction is strong at lower dimensions, sev- 
eral recent experimental investigations on the spin state 
of OD H and ID systems ^ have indicated the pos- 
sibility of a SSP, even though the origin and nature of 
such a phase is controversial. In low-density 2D elec- 
tron systems (2DES's), however, a direct observation of 
a "spin gap" has remained experimentally illusive, even 
though evidence of an enhanced g-factor and anomalous 
spin-susceptibility have been reported when the disorder 
is low ,7,, 1^ ^ ilOiJ . 

It is now known that, an unpaired spin in weakly- 
isolated systems are screened by Kondo-like singlet for- 
mation with the electrons in the leads Hdll. This re- 
sults in a resonance in the tunnelling density-of-states 
(DOS), and an enhancement in the differential conduc- 
tance (dl/dV) through the system when the chemical po- 
tential of the leads are aligned. When studied as a func- 
tion of source-drain bias {Vsd), the zero-bias peak (ZBP) 
splits linearly in finite B by Az — g* ^bB, where g* and 
/is are the effective g-factor and the Bohr magneton re- 
spectively. Recent non-linear studies in ballistic quantum 
point contacts (QPC's) have also shown the evidence of a 
ZBP J^], which was attributed to a Kondo-like correla- 
tion resulting from a dynamic SSP in ID. In this Letter, 
we report the direct observation of SSP in low-density 
2DES's. High quality 2DES's of mesoscopic length scales 
were used with no intentional in-plane confinements. In 
the localized regime, clear evidence of a split ZBP was ob- 



served at B = 0. The magnitude of the split (A) evolved 
continuously with B, implying its origin to be related to 
the underlying spin structure of the 2DES. Simultaneous 
measurement of the Fermi energy {Ep) showed that such 
a polarization is an intrinsic ground state property of 2D 
systems, and depends critically on the impurity scatter- 
ing, as well as the temperature (T) and magnetic field 
{B). 

We have used 2DES's confined to the triangular poten- 
tial wells of Si (5-doped GaAs/AlGaAs heterostructures. 
In most samples, the 2DES's were formed « 300 nm be- 
low the surface. Disorder was varied by changing the 
spacer thickness Sgp separating the dopant layer from the 
GaAs/AlGaAs interface. Data from two samples with Sgp 
= 40 nm and 60 nm (referred as A78 and A79, respec- 
tively) are reported in this work, even though all the sam- 
ples show qualitatively similar results. In both samples, 
a source-drain voltage Vsd was applied on a 5/j,mx5/j,m 
region of the wafer, defined by an etched mesa and a 
metalhc surface gate. At zero gate voltage {Vg = 0) 
the mobility of both samples were ^ 2 x 10^ cm^/V- 
sec. Lifetime measurements confirm the dominance of 
the small angle scattering. By varying Vg, electron den- 
sity {ris) as low as ~ 5 X 10^ cm~^ could be attained 
in sample A79 (corresponding to an interaction param- 
eter Ts = l/a*Qy/Tms 7.6, where a*g is the effective 
Bohr radius). In all the magnetic field measurements, 
B was applied in the plane of the 2DES and parallel 
to the direction of the current. The differential con- 
ductance dl/dV was measured with a standard 2-probe 
mixed ac-dc method, where the ac excitation bias was 
kept at < ksT/e. For all samples, the gate was trained 
several hundred times to obtain an excellent run to run 
reproducibility (better than 0.1%). 

Dependence of the linear response conductance (G) on 
Vg is shown in Fig. 1 at various values of B, recorded 
at T « 35 mK in sample A79. We focus on the local- 
ized regime (MIT in A79 occurs at G ~ 3 x e^/h), where 
the electron density is low, and the interaction effects are 
most pronounced. For all VJ^, G was found to decrease 
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rapidly with increasing B. The functional dependence of 
the in plane magneto-conductance (MC) on B has been 
shown to change from ^ e^^ to ^ e^^ at a critical field 
Be, when the electron gas becomes completely spin polar- 
ized, and the majority and minority spin bands are sepa- 
rated by the Fermi energy ^ J^Qj . In the inset of Fig. 1 we 
have shown the zero-bias MC as a function of B^ at four 
representative gate voltages. The values of Be at each Vg, 
shown by the arrows, indicate a deviation from the lin- 
ear behavior at low fields. This enables us to evaluate the 
spin non-degenerate Fermi energy Ep as Ep = g*^BBc- 
The dependence of Ep on Vg is shown in Fig. 2c. We have 
used a 5* = 3.4 1 (?h I for evaluating Ep, where \gh\ = 0.44 
is the band g-factor in bulk GaAs. This value of g* is 
obtained from Vsd measurements, and will be discussed 
later. Note that above Vg « —0.376 V, corresponding 
to Tis ~ 5 X 10® cm~^, the transport is essentially 2D 
in nature with E*p varying approximately linearly with 
Vg. The slope dE*p/dVg was found to be « 12 meV/V, 
agreeing roughly with the effective free-electron spin 
non-degenerate 2D density-of-states, {dE*p / dVg) free = 
{h'^/2TTm*)dns/dVg « eoerh"^ /2nem*ds ~ 16 meV/V, 
where dg = 310 nm is the depth of the 2DES from the 
surface, and m* w 0.067me is the band effective mass of 
the electron. The discrepancy could be due to a weak 
density-dependence of g* 0- Below Vg « -0.376 V 
(G « 0.3 — 0.4 X e'^/h), Ep drops abruptly, possibly due 
to the onset of inhomogeneity in the charge distribution 
as screening becomes weak. 

The Vsd dependence of dl/dV at various values of 
Vg (i.e.. Us) is shown in Fig. 2a. Vg differs by 1 mV in 
successive offset-corrected traces. The striking feature of 
these traces is the double peak structure of dl/dV with 
a local minimum at Vsd — 0. This was found to be a 
generic feature observed in all the low-disorder samples 
of similar dimensions. The detection of the double-peak 
structure was difficult in the strongly localized regime 
(Vg < —0.385 V), as well as in the metallic regime 
{Vg > —0.35 V), indicating a disorder-dependent window 
of Us where the effect becomes clearly visible. In most 
cases the peaks are dissimilar in magnitude and width, 
both of which vary when Vg is changed. The separation 
(A) of the peaks shows a non-monotonic dependence on 

Vg, as shown in Fig. 2b. At low Vg ( 0.385 V), A is 

largest, but decreases rapidly with increasing Vg, reach- 
ing a minimum at Vg ~ —0.375 V. Comparing to the Vg 
dependence of Ep (Fig. 2c) we find the onset of linear de- 
pendence of E}. at the same Vg. When Vg > -0.375 V, A 
increases roughly linearly with increasing Vg. Since this 
regime can be directly associated with a 2D ground state, 
we shall restrict further discussions on A to this range of 
Vg. Extrapolating the linear dependence of A, we find 
the Vg (« —0.395 V) at which A = 0, agrees within 
the experimental uncertainty with the Vg (« —0.408 V) 
at which the extrapolated Ep = 0, establishing a direct 
correspondence between A and n^. 



In order to investigate the nature of this effect, we 
have then studied the double peak structure in a parallel 
magnetic field. The behavior is illustrated with the trace 
measured at Vg « —0.37 V, indicated by the arrow in 
Fig. 2b. This is shown in Fig. 3, where the traces with 
increasing B are vertically offset for clarity. The peak 
positions (V^) change non-monotonically as a function of 
B (see the arrows). We find that the peaks close in over 
the field scale oi B < 0.5 T, but eventually separate at 
higher fields (see the inset of Fig. 3). The evolution of 
the peaks with B strongly suggests a spin-related effect, 
where A is the energy difference between the opposite 
spins. Note that the i3-field does not split each peak 
individually, differentiating our case from a disorder in- 
duced quantum molecule. An asymmetrical suppression 
of the peaks with increasing B sets the maximum field 
scale (« 2.8 T) of our experiments, beyond which the left 
peak becomes undetectable. 

The observability and characteristics of the peaks were 
found to depend critically on disorder and T. We have 
changed the local disorder profile over the sample region 
by controlled thermal cycles from room temperature to 
4.2 K. The result for three successive thermal cycles car- 
ried out for A78 are shown in Fig. 4a. All traces (B and 
C are shifted vertically for clarity) were recorded at a 
similar Ep and B = 0. While the general behavior of 
dl/dV is similar and agrees with that of A79, we find 
the overall magnitude and width of the peaks to differ 
markedly, even from one cooling to the other. In gen- 
eral, when the cooling was done at a slow rate (trace A: 
over a few hours) we found the peaks to be more pro- 
nounced than when the cooling was done rapidly (trace 
C: over few tens of minutes). Greater disorder in C is 
also observable in terms of the linear response conduc- 
tance, G « O.le^/h, which, in case of A is ~ 0.3e^//i. 
Disorder broadening also affects the observability of the 
effect as the sample dimensions are increased. Typically, 
no double-peaks were observed in 10/zmxlO/zm samples. 

Thermal broadening of the peaks, associated with a 
strong suppression of the peak height, from T ^ 70 mK 
to ~ 0.5 K are shown in Fig. 4b. Note that even though 
G a.t Vsd — rises with increasing T, as expected in a 
localized system, it is essentially a result of the overlap 
of the broadening peaks. This also confirms that as a 
function of Vsd, we indeed observe a non-linear effect, 
and the double peak structure is a ground-state property 
of the 2D electron system in the low density regime. 

We now discuss possible mechanisms which may give 
rise to such a structure in non-equilibrium measurements. 
Assuming the formation of a weak tunnel barrier, which 
physically splits the 2D region in two parts, the conser- 
vation of the transverse momentum will allow tunnelling 
through the barrier only at Vsd = Q- This would result 
in a disorder-broadened ZBP in the tunnelling dl/dV. 
An interaction- induced suppression of states at Ep, e.g., 
~ In (|e|r) in the diffusive regime 13] or Efros-Shklovskii- 
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type soft gap |e| in the hopping regime may 
then spil t the ZBP to give rise to the double peak struc- 
ture [l^- The main arguments against such a scenario 
are, (a) the separation of the peaks are -B-independent, 
hence inconsistent with Fig. 3, and (b) the separation 
decreases with increasing Ep, contrary to the result of 

Fig. 2b Hi. 

To investigate if the peak separation A is indeed spin- 
related, we have plotted A as a function ofB for the data 
shown in fig. 3. As shown in Fig. 5, at low fields A de- 
creases with increasing field, while at high i? (> 1 T), A 
shows monotonic increase as B increases. In this high-_B 
regime, A tends asymptotically to a linear B dependence 
(the solid line), which when extrapolated, passes through 
the origin. This we identify as Zeeman spin splitting. 
From the asymptote, Az = g*/JBB, we find g* to be 
« 1.49 = 3.4|(7f,|, agreeing closely to the values reported 
in recent measurements over similar range of H. The 
peak positions thus represent the energies of two oppo- 
site spin states in a magnetic field, and from the non-zero 
separation at _B = 0, signify a non-zero spin polarization. 

From the critical role of spin , and also the behav- 
ior of dl/dV near Vsd = 0, a Kondo-like many-body 
correlation, as discussed extensively in the context of 
quantum dots can be envisaged. There are how- 
ever important differences. In our system, as in the case 
of QPC's, there are no obvious singly-occupied quasi- 
bound electronic states. Unlike the suggested ferromag- 
netic states in QPC's the magnetic moment of a 
frozen spin-polarized state in 2D would be much higher 
than that of a single electron, and hence Kondo-type 
screening would be difficult. However, as indicated in 
Fig. 2a for Vg > —0.375 V, the spin-polarization in our 
system seems to be dynamic, and increases with Ug- This, 
as suggested for QPC's 0, in presence of optimal band- 
hybridization, could result in a Kondo-enhanced ZBP in 
dl/dV. Our case is however further complicated by the 
splitting of the ZBP at _B = 0, resembl ing the behavior of 
ZBP in dots and QPC's at a finite B |ll[i3|. Although 
the Kondo effect in coupled dot systems has a similar 
non-equilibrium behavior as a function of Vsd at _B = 0, 
such a case would appear only over a restricted parameter 
range of inter-dot and lead couplings, and hence should 
be rather rare in open systems with no intentional con- 
finement 10. 

Finally, in parallel B, it has been shown that in suffi- 
ciently smooth disorder, both SO and exchange-induced 
spin splitting could result in satellite-peaks in the tun- 
nelling DOS at finite bias w ±y/A^ + A|^, where Aq is 
the magnitude of SSP |ia]. The tunneUing DOS could 
be obtained in a non-linear dl/dV measurement if we as- 
sume a quasi-ballistic transport in our clean mesoscopic 
samples, i.e., electron energies are lost mainly at the 
leads. Apart from explaining the double-peak structure 
in dl/dV, this also justifies, (a) broadening of the peaks 
with disorder and (b) the linear i?-dependence of A at 



high B, where A —^ Az- If Ag arises from a strong SO 
coupling, one expects negative magneto-conductance at 
low perpendicular iJ-field, arising from weak antilocal- 
ization. No evidence of antilocalization was observed in 
our mesoscopic 2D systems. The SO-origin of Aq seems 
to be unlikely on two more grounds. First, calculation of 
the magnitude of splitting in GaAs/AlGaAs heterostruc- 
tures, taking into account both bulk and Rashba terms, 
shows that Aq would be ^ 0.05 meV below Ug ^ 10^^ 
cm~^ Secondly, if the absence of antilocalization is 
attributed to a B-dependent Aq, as in large open quan- 
tum dots, Aq increases with B at low fields |a,|2l|. How- 
ever, as shown in the inset of Fig. 5, Aq ~ a/A^ — A^ 
decreases with increasing B in our samples. Alterna- 
tively, if we assume that Aq arises due to exchange, and 
in the Fermi liquid limit, represents the shift of one spin 
state with respect to the other (see the inset of Fig. 5), 
then an estimate of polarization C could be obtained as 
C — (rtsj — nsi)/ns ^ A/2Ep. In the regime of 2D con- 
duction, Vg > —0.375 V, we find ( ~ 0.1, implying only a 
partial spin polarization. This is probably not surprising 
considering the low Ts ~ 6 — 7 in our system, and finite 
T Q. The exchange-origin of Aq is also supported by 
the decrease of Aq with increasing J5, as the additional 
confinement (^ B^) imposed by B tends to counteract 
the formation of parallel spins [la| ■ 

In conclusion, non-linear conductance measurements 
in high quality mesoscopic 2D electron systems show di- 
rect evidence of spontaneous spin polarization at low 
electron densities. This spin polarization is density- 
dependent, and manifested by a split zero-bias peak 
in source-drain bias at i? = 0. The peaks are gener- 
ally asymmetric and strongly depend on the local dis- 
order. Simultaneous measurement of the Fermi energy, 
and temperature-dependence of conductance confirm the 
splitting to be an intrinsic property of 2D systems, and 
also identify a window of electron density and temper- 
ature over which the splitting becomes experimentally 
observable. 
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FIG. 1: Typical gate voltage Vg dependence of the linear response conductance G at various magnetic fields. Inset: G, as a 
function of at four V^'s. Arrows indicate the deviation from linear dependence, and the onset of complete polarization. 



Fig. 2 




FIG. 2: (a) The differential conductance dl/dV vs. source-drain bias Vsd at various Vg. Vg differs by 1 mV for successive 
sweeps, (b) The peak separation A as a function of Vg. (c) Vg-dependence of the spin non-degenerate Fermi energy _EJ, 
obtained from the magnetic field scale of complete field polarization (see text). 
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FIG. 3: Evolution of the differential conductance with magnetic field B. We have illustrated this with the trace obtained at 
Vg w —0.37 V (indicated by the arrow in Fig. 2b). Traces are shifted vertically for clarity. Inset: Positions (Vp) of the peaks 
as a function of B. 



Fig. 4 
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FIG. 4: (a) Double-peak structure in dl/dV in the same sample for various cool-downs. Traces arc shifted vertically for easy 
comparison. The linear response conductance is minimum for C and maximum for A. (b) Temperature dependence of the 
double peak structure. Note the smearing of each peak results in the increase of the zero-bias conductance with increasing T. 
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FIG. 5: Dependence of the peak separation A on the magnetic field B. We have used the A shown in Fig. 3. The high-field 
asymptote extrapolates through the origin. Inset: B-dependence of the Ao. 



